Charge Friedel oscillations in a Mott insulator 
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When a metal undergoes a transition to an insulator it will lose its electronic Fermi surface. 
Interestingly in some situations a 'ghost' Fermi surface of electrically neutral spin carrying fermions 
may survive into the insulator. Such a novel ghost Fermi surface has been proposed to underlie 
the properties of a few different materials but its direct detection has proven elusive. In this paper 
we show that the ghost Fermi surface leads to slowly decaying spatial oscillations of the electron 
density near impurities or other defects. These and related oscillations stem from the sharpness of 
the ghost Fermi surface and are direct analogs of the familiar Friedel oscillations in metals. The 
oscillation period contains geometric information about the shape of the ghost Fermi surface which 
can be potentially exploited to detect its existence. 



I. INTRODUCTION 

In recent years it has become clear that the insulating 
side in the vicinity of the Mott metal-insulator transition 
may, in some frustrated lattices, provide a realization of 
the long sought quantum spin liquid state. These are 
insulators with an odd number of electrons per unit cell 
and that do not order magnetically or break lattice sym- 
metries which would double the unit cell. Examples are 
the layered quasi-two dimensional organic materials 1,2 
k-(ET) 2 Cu 2 (CN) 3 and EtMe 3 Sb[Pd{dmit) 2 ] 2 , and the 
three dimensional 'hyperkagome' material^ Nailr 3 0%. 
These materials can all be driven metallic by applica- 
tion of moderate pressure^, and hence are appropri- 
ately thought of as weak Mott insulators which are not 
too deep into the Mott side of the metal-insulator tran- 
sition. Indeed quantum spin liquid behavior may well 
be a common fate of weak Mott insulators. In contrast, 
the quantum spin liquid state remains much more elu- 
sive in strong Mott insulators with geometric frustration. 
A notable exception is the spin- 1/2 Kagome magnet® 
ZnCu 3 {OH) & Cl 2 . 

A fundamental question about the Mott transition is 
the fate of the Fermi surface of the metal as it transitions 
into the Mott insulator. Theoretically the natural spin 
liquid state in the vicinity of the Mott transition retains 
a 'ghost Fermi surface' of mobile neutral spin carrying 
fermionic excitations (dubbed spinons) even as the charge 
of the electron gets pinned to the lattice to form an insu- 
lator. Such a state has been proposed to describe the in- 
termediate temperature physics of many of the candidate 
spin liquid materials 7 ^. Experimentally both organic 
salts and Na^Ir^Os have gapless spin excitations with a 
low temperature specific heat^ and spin susceptibility 1,2 
resembling that of a metal even though they are insulat- 
ing for electrical transport. Further recent experiments 11 
in EtMe 3 Sb[Pd(dmit)2\2 show metallic thermal trans- 
port in this insulator, suggesting that the gapless excita- 
tions are, in fact, mobile. 

How can we really tell if such a ghost Fermi surface of 
spinons survives in a weak Mott insulator? If the spinon 
Fermi surface state survives to low-T, Motrunicbi 2 - has 
discussed the possibility of unusual quantum oscillation 



phenomena in a magnetic field to detect the Fermi sur- 
face. However in practice it is likely that the spinon 
Fermi surface undergoes an instability at low-T, possibly 
related to spinon pairing^—. Consequently other less 
direct methods to detect it at intermediate T are called 
for. Another idea 16 is to look for an oscillatory coupling 
between two ferromagnets separated by a spin liquid ma- 
terial due to RKKY spin correlations associated with the 
spinon Fermi surface. However unlike the standard giant 
magnetoresistance effect the oscillatory coupling cannot 
be readily detected by measuring the electrical resistance, 
and some further idea is required. This scheme also re- 
quires controlling the thickness of the spin liquid layer to 
atomic precision, and hence is challenging. 

In this paper we identify new and surprising aspects 
of such a spin liquid Mott insulator that may potentially 
offer a simpler route to detecting the spinon Fermi sur- 
face. First we show that the electrical charge density in 
response to an external potential has Friedel oscillations 
at the "2_Kf" wavevectors that connect tangential an- 
tipodal portions of the spinon Fermi surface. Friedel os- 
cillations in the charge density are a well known property 
of metallic systems with a sharp electron Fermi surface. 
Remarkably these oscillations survive the transition into 
the spin liquid Mott insulating state. The spinons that 
form the Fermi surface in this spin liquid carry spin, and 
hence there are 2K-p singularities in the spin response. 
However spinons are electrically neutral - despite that the 
2Kp singularities also show up in the charge response. 

Fundamentally the charge Friedel oscillations stem 
from 2K-p singularities in the spinon number density re- 
sponse. This is a spin singlet response and should be dis- 
tinguished from the spin density response which is a spin 
triplet. In a weak Mott insulator, short wavelength fluc- 
tuations of the spinon density will couple to the electrical 
charge density and lead to the charge Friedel oscillations. 

A closely related effect is the existence of a Kohn 
anomaly at the 2K-p wavevectors in the phonon spectrum 
arising from the coupling of the phonons to the spinon 
density. Another related effect will appear in Scanning 
Tunneling Microscopy (STM) spectra above the thresh- 
old bias for tunneling into the Mott insulator. Near the 
threshold the tunneling spectrum will show spatial mod- 



ulation near impurities with wavelength set by the 2K-p 
wavevectors of the spinon Fermi surface. These related 
effects are more directly measurable in experiments and 
can serve to detect the spinon Fermi surface. 

We provide an estimate of the magnitude of the oscil- 
lation amplitude of the charge density and show that it 
is essentially unchanged as the material undergoes the 
metal-insulator transition. However it is rather weak 
compared to the free Fermi gas so that its direct mea- 
surement may be even more challenging that in a nor- 
mal metal. Not surprisingly, the oscillation magnitude 
decreases with increasing Mott gap, eventually going to 
zero in the extreme limit of a spin model with no charge 
degrees of freedom. In contrast the Kohn anomaly is 
not similarly weakened and survives even into the strong 
Mott insulator limit. Just as in a normal metal, it is 
therefore likely easier to detect the Friedel oscillations 
through their effects on the phonon and STM spectra. 
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II. SPINON FERMI SURFACES AND 2K F 

SINGULARITIES IN A WEAK MOTT 

INSULATOR 



For concreteness consider a one band Hubbard model 
at half-filling on a non-bipartite lattice such as the trian- 
gular lattice: 



H = -t J2 (4 a Cr> a + h.c)+Uj2(nr-l) 2 (1) 

<rr'>a r 



where c ra destroys an electron with spin a = J [, X at site 
r. n r = c\c r is the electron number at site r. U > is 
an on-site repulsion. As the ratio g = t/U decreases the 
system undergoes a Mott transition from a Fermi liquid 
metal to a Mott insulator. We are interested in situa- 
tions in which this Mott insulator is a non-magnetic spin 
liquid with gapless spin excitations. The Mott transition 
and the spin liquid phase are conveniently discussed us- 
ing the slave rotor representation of Ref. [lj. We write 
c ra = e'^/ ra with e l ^ r = b r a spin-0 charge-e boson, 
and f ra a spin-1/2 charge-0 fermionic spinon. The elec- 
trical charge at a site n r is identified with the boson 
number and is conjugate to the boson phase <p r . The 
physical electron operator is manifestly invariant under 
local opposite phase rotations of b r and f ra respectively. 
A proper reformulation of the Hubbard model in terms of 
the b r and f ra will thus include an emergent U(l) gauge 
field. 

An effective theory 8 for the phases of the Hubbard 
model is given by an imaginary time path integral with 



Cb = ^2 -Tj- 1 (d T + ia a + A) b r \ 2 - ^2 tb (b*b' r e arr ' + < 
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(3) 



£/ =^2fr(d T - (ia + X))f r - ^2 tf(frfr'e ' a -' + c.cj 

r <rr'> 

(4) 

Here ao,a rr i are the temporal and spatial components 
of the emergent U(l) gauge field. The parameters %j 
are boson and spinon hopping amplitudes which may be 
estimated in a mean field approximation. A is adjusted 
to ensure that (rib — n/} = 0. 

In a mean field approximation, which ignores the gauge 
but not other interactions, the charge sector is mapped 
to a boson Hubbard model with hopping tb and boson 
repulsion U while the spin sector is described by free 
fermionic spinons that form a Fermi surface. When t/U 
is large, the boson hopping tb dominates leading to a state 
with (b r ) ^ 0. In terms of the original electrons this is 
a Fermi liquid phase. As t/U is decreased, the bosons 
form a Mott insulating state with a charge gap. This 
corresponds to a spin liquid electronic Mott insulating 
state which retains a Fermi surface of gapless neutral 
spin-1/2 spinons. Beyond mean field, at low energies in 
the Mott insulating phase the bosons may be integrated 
out, resulting in an effective theory of a spinon Fermi 
surface coupled to a U(l) gauge field. 

The possibility of sharp 2Kf singularities for spinons 
was addressed in an early worbi£. Very recently con- 
trolled calculations 19 of 2Kp singularities in both the 
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spinon number rlf r = fjf r and the spin S r 
been performed. Consider for instance the static density 
response function Xf (p) °f the spinons at a wavevector 
p. If Q is any 2Kp wavevector that connects two antipo- 
dal tangential portions of the spinon Fermi surface then 
X/(Q + q) is singular for small |q|: 

Xj )Q (q) = X /(Q + q) - X/(Q) = k|| 1^ (j^f) , (5) 

where gn , q± are the magnitudes of the components of 
q that are parallel and perpendicular to Q respectively. 
The exponent <f> and the scaling function F depend 
on spatial dimensionality. For instance in two dimen- 
sions and in the absence of the gauge field the exponent 
(f) = 1/2. The presence of the gauge field modifies the 
exponent in a manner that is described in Ref. Il9l . In 
three dimensions the gauge field plays a more innocuous 
role and the exponent is the same as in the usual three 
dimensional Fermi liquid. Thus the spinon density cor- 
relations have sharp singular structure at all such 2A"p 
wavevectors of the underlying spinon Fermi surface. Ex- 
actly the same exponent <f> and scaling function F also 
describe the 2K-p singularities in the spin density correla- 
tions due to an emergent symmetry of the spinon Fermi 
surface stated. 



III. 



CHARGE DENSITY CORRELATIONS 



We now show that in a weak Mott insulator the univer- 
sal short distance 2Kp oscillation in nj r imprints itself 
on the correlations of the physical electrical charge n r as 
well. Further the magnitude of the oscillations are more 
or less unchanged across the metal-insulator transition. 

The spinon number nf >r is a spin singlet operator that 
transforms under lattice and time reversal symmetries in 
the same way as the electrical charge n r = c\c r . Thus on 
general symmetry grounds we might expect that there is 
a linear coupling between nj and n so that any structure 
in the nf correlator is transferred to n. For long wave- 
length external perturbations that couples to the charge 
density (which is the boson density), the boson system 
is gapped and the charge response is that of an incom- 
pressible insulator. However for short wavelength pertur- 
bations the boson density will have a non-zero response 
which, in turn, induces a response of the spinon den- 
sity. If the wavevector matches a 2K? wavevector of the 
spinon Fermi surface, there will be long range Friedel os- 
cillations of the fermion density that couple back to the 
physical charge density which will therefore show Friedel 
oscillations. 

An approximate self-consistent calculation of the mag- 
nitude of this effect may be done within the slave parti- 
cle mean field theory. An external potential with Fourier 
components V(p) induces some change SX in the chem- 
ical potential. In general the hopping parameters £;,,£/ 
will also change, for a rough estimate we ignore these 
below. The resulting change in the boson and fermion 
densities may be expressed in linear response theory as 



S(n)(p) 
S(n f )(p) 



Xb(p)(V(p) + 
-X/(p)<5A(p) 



<JA(p)) 



(6) 

(7) 



Here Xbj are the density response function of the bo- 
son and fermion respectively. The condition S(n)(p) = 
<5(n/)(p) determines SX. The charge density response 
function defined through 5(n)(p) — x(p)^(p) then takes 
the form 



X 



XbXf 
Xb + Xf 



(8) 



This is the well-known Ioffe-Larkin composition rule 
which we now use to evaluate the 2Kp response. 

Near a 2Kp wavevector p rs Q, the spinon density 
response Xf has the singularity discussed above. How- 
ever as we move across the metal-insulator transition, 
though the small q response of the boson change strik- 
ingly the high q response will be only weakly modified. 
We write the fermion density response X/(Q + Q) = 
Xf ,2K F + Xf ,2K F (q) where Xf , 2 k f is the smooth back- 
ground response at 2Kp and xf 2K (l) ^ s the singular 
part. Denoting x&(Q + q) = X&: we get for the singular 
part of the electrical density response 



where the reduction factor R is given by 
R = 



1 



l /2K F 

"TO — 



(10) 



X2K F 



Rx s f , 



2Kb 



(9) 



R is the amount by which the Friedel oscillations of the 
electrical charge density are reduced compared to the 
oscillations of the spinons. As jj is varied across the 
Mott transition both Xb an d X ? ik y wm " be essentially 
unchanged so that the factor i?, and hence the magni- 
tude of the charge Friedel oscillations, is more or less the 
same on either side of the transition. R may be estimated 
within the slave rotor mean field theory coupled with a 
large- A approximation to the quantum rotor model that 
describes the bosons (see Appendix |A"| . Near the Mott 
transition, we find R = 0.024. Though the magnitude 
of the charge Friedel oscillations is thus very small it is 
interesting conceptually that it is there at all. In the 
next section we will see how we may directly probe in 
experiments the 2K^ singularities of the spinon density 
response. Going deeper into the Mott insulator decreases 
R . The extreme limit of a spin model corresponds to 
Xb(q) = for all q and R = as expected. 

Actually as discussed above the detailed nature of the 
2K~y singularity for the spinons is modified and even po- 
tentially enhanced^ in the Mott insulator as compared 
to the Fermi liquid due to the coupling to a gapless U(l) 
gauge field. Thus we reach the remarkable conclusion 
that not only are there charge Friedel oscillations in the 
insulator but they may even be more slowly decaying 
than in the metal. 



IV. DETECTION OF SPINON DENSITY 
FRIEDEL OSCILLATIONS 

The presence of the 2Kp singularities in the charge 
density response relies fundamentally on the 2Ap singu- 
larities in the spinon density response. We now argue 
that this latter singularity directly affects the phonon 
and STM spectra which could be more directly useful to 
detect the spinon Fermi surface. 

In a metal the singular 2Kp density fluctuations im- 
print themselves on the phonon spectrum through the fa- 
mous Kohn anomaly. Here we show that a similar Kohn 
anomaly occurs in the spin liquid Mott insulator. In the 
spin liquid, the charge of the electron gets pinned to the 
ions that make up the underlying lattice, while the spin 
continues to be mobile and forms a Fermi surface. Thus 
the oscillating entity in a phonon mode is the electrically 
neutral ion-boson composite (see Figs. Ill2p . However 
this neutral object carries gauge charge of the boson b r , 
and hence there will be long range interactions between 
different ion-boson composites. Consequently to prop- 
erly describe the phonon spectrum we must take into ac- 
count the screening of the emergent gauge interaction by 
the spinon fluid. Specifically there will be an emergent 
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FIG. 1: In a regular metal, charge (represented by white half- 
spheres) and spin (black half-spheres with arrows) of the elec- 
trons are confined, and the electrons move to screen any dis- 
placed positively charge ions (big spheres). 



FIG. 2: In the spin liquid Mott insulator, the charge of the 
electron is bound to the ions while the spins stay mobile. 
The electrically neutral ion-chargon composite carry internal 
gauge charge however and their displacement is screened by 
the (oppositely) gauge-charged spinon fluid. 



'Coulomb' interaction between the ion-boson composites 
that will be screened by the mobile spinons. As in the 
usual discussions of the Kohn anomaly the ability of the 
spinon fluid to screen the gauge interaction is sharply 
changed as the phonon wavevector passes through any 
2Kp wavevector. Consequently there is a Kohn anomaly. 
This argument also makes clear that the Kohn anomaly 
will survive even deep into the Mott insulator so long as 
the system retains the spinon Fermi surface. Thus it is 
not weakened by the suppression of electric charge fluctu- 
ations, and will be more or less unchanged from a weakly 
interacting metal. 

In metals Friedel-like oscillations have been imaged 
through Scanning Tunneling Microscopy. The most com- 
mon strategy is to measure the spatial modulation of 
the tunneling spectrum at fixed bias. This modulation 
comes from the standing wave pattern caused by interfer- 
ence between plane wave states that are mixed by defects. 
Strictly speaking these are distinct from the Friedel oscil- 
latons but at low bias the modulation is again at the 2Kp 
wavevectors of the Fermi surface. In the Mott insulator 
we are faced with the immediate problem that (at T = 0) 
there is a gap in the single electron spectrum. However 
if the bias exceeds the gap it will be possible to tunnel 
into the sample. We show below that the tunneling spec- 
trum at a bias voltage close to the gap will acquire spatial 
modulation at the 2Kp wavevectors of the spinon Fermi 
surface. The tunneling spectrum near the threshold is 
readily calculated within the slave rotor mean field the- 
ory. The local tunneling density of states N(uj) is then a 
convolution of the boson and fermion densities of states 



N(u) = / d£lN b (fi) N f (u - n) (9((Q) - (Q - u)) 

(11) 

In Appendix IB1 we show that N(ui) = for \ui\ < A, and 
increases linearly just beyond the threshold, i.e 



TABLE I: Smallest wave-vectors connecting antipodal tan- 
gential portions of the Fermi surface in the spin-liquid candi- 
dates 



K - {ET) 2 Cu2{CN)3 



Triangular 
0.58- « 0.22 A 



ZnCu 3 (OH) 6 Ck 



Kagome 
0.58- « 0.265 A 



NailrsOs 



Hyperkagome 
hole: 0.52^ « 0.18 A 

a 

el.: 0.62^ « 0.21 A 



iV( w ) = ^(M-A)|(M-A|). 



(12) 



In the presence of a defect such as a step edge or an 
impurity, Nf(u>; x) oscillates as a function of spatial loca- 
tion x at the 2Kp wavevectors of the spinon Fermi surface 
for low frequency u — >• . Nb(w) will also have some x 
dependence but will not oscillate at the same wavevec- 
tors. Near the threshold the tunneling electron creates a 
boson with energy close to A and a fermion close to zero 
energy. As shown in Appendix [B] the coefficient A will 
show spatial modulation at the 2Kp wavevectors. 

We now comment on the observability of these 
effects in the candidate spin materials in hand, 
namely the quasi-two dimensional k — {ET)2Cu2{CN)z 
and EtMe^Sb[Pd{dmit)2]%, and the three dimensional 
Na^Ir^Os- In all three materials determination of 
the phonon spectrum (to look for the Kohn anomaly) 
through neutron scattering is challenging but may pos- 
sibly be done through inelastic X-ray scattering tech- 
niques. In the organic materials the spinon Fermi tem- 
perature is estimated to be ~ 250A'. Consequently the 
Kohn anomaly will only be visible at much lower temper- 
ature. Furthermore due to the layered structure it will be 
independent of the component of momentum perpendic- 
ular to the layers. Some representative wavevectors along 
high symmetry directions are given in Table HI and the 
surface of all 2Ap wavevectors in the first Brillouin zone 
is shown in Fig. IIVI assuming a perfect triangular lattice. 
The hyper- kagome material Na^Ir^Os has a rather com- 
plicated band structure in its metallic state which will 
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Appendix A: Mean field theory in the large N limit 

In the mean field theory, the fermion action is gaus- 
sian but the boson action is still interacting. We gen- 
eralize the 0(2) model a model of N complex bosons 
b = (bi, . . . , 6jv) which we solve by expanding in pow- 
ers of 1/N and extrapolating to N = 1. The mean-field 
boson action is 



Cb = zZ^j U + Y,Ay i + Ahi 



(Al) 



FIG. 3: 2kF vectors for the triangular lattice at half filling in 
the 1st Brillouin zone. 



J2 *6 (K ■ b r' + c - c ) + iK(\b r \ 2 - A), (A2) 



be inherited by the spinons in the proposed spin liquid 
Mott insulator. Some representative 2 Ky wavevectors 
were calculated for this material in Ref. [la, and we list 
them here for completeness in Table |U 

Previous STM studies^ of the n — ET organics have 
focused on a superconducting material. STM studies of 
the insulating spin liquid will be useful both in determin- 
ing the single particle gap and in detecting the possible 
spinon Fermi surface. In this context we note that in 
K — (ET) 2 Cu2(CN) 3 the charge gap as measured by op- 
tical conductivity and dc transport is small (in the 15 
-100 meV range) though it is hard to extract a precise 
number. This range should be readily accessible in STM 
experiments. 

Finally we comment briefly on strong Mott insula- 
tors with negligible virtual charge fluctuations but which 
nevertheless are in quantum spin liquid phases. This 
is the situation, for instance, in the Kagome material 
ZnCu^(OH)eCl2- Such systems are well described by 
simple spin-only models with dominant nearest neighbor 
exchange. Spin liquids with spinon Fermi surfaces may 
emerge in these spin-only models as well. Indeed such 
a spin liquid has been proposed for ZnCu^(OH)(,Cl2 in 
Ref. [2l|. What is the interpretation of the 2K-p singular- 
ities in the spinon number n/ in that case? As this is a 
spin singlet this singular 2Kp response signals slow power 
correlations of incommensurate valence bond solid order. 
Such correlations have been found in a spin-1/2 system 
on a triangular strip22. The arguments presented above 
show that phonons of the lattice will couple to these in- 
commensurate valence bond fluctuations (similarly to the 
familiar one dimensional gapless spin liquids), and there 
will again be a Kohn anomaly in the phonon spectrum 
at the 2Kp wavevectors. Detecting such an anomaly will 
then be a non-trivial proof of the proposal that a spinon 
Fermi surface exists in this material. 



where g l are the 2A 2 + 0(A) generators of 0(2N) and 
A is a Lagrange multiplier field, that is integrated over in 
the path integral. The external gauge field A l should not 
be confused with the emergent gauge field a^ in Eq. ([2")l, 

A 2 

For large N we can replace ^j = iX by its saddle-point 
value, which for g < g c is given by 



duj 



Z^ 

fceBZ 



I 9 A2 1 



(A3) 



2C 



2/ 



where e^ = 6 — ^ e M i fc ; with Si all vectors connecting 
a lattice site to its nearest neighbors. At half filling on 
the triangular lattice tb ~ 0.33i and by setting A = we 
find g c fa 0.39. The boson-susceptibility \b = J2i xl = 
X& la + X& ara where the diamagnetic and paramagnetic 
contributions are given by 



para / -*. 

Xb (q) 



dia 2N 2 



2N 2 
(Ut b ) 



dw Yl 



, ^2 + A 2 C^ +A 2 

fcesz k ^ 2Ut b fc fc+g^ 2Ut h 



-xr a (o). 



(A4) 



The fermion susceptibility is given by the usual expres- 
sion 



XI 



(q) = N s JdCo J2 



1 



1 



fceBZ 



-iu + t/eg -n-iu + t f e % g - n 



(A5) 

with N s the number of spin species. For N = 1, N s = 2, 
A = 0, and q = q Q x, we find Xb(qo) ~ ^)Xf(?o) ~ TT^ 
so that 



\ Xb(qo) 



0.024 



(A6) 



Appendix B: Tunneling d.o.s. 

The electron Green function in real space is given by 

G e (T,r) = G f (r,r)G b (r,r), (Bl) 

so that 



1 



N{u) = Im / d 2 kG r e et (uj,k) 



(B2) 



dnN b (O) N f (u - n) (9 (O) - (fi - w)) , 

(B3) 



where the Heaviside-functions 9 are a consequence of the 
pole structure of the retarded Green functions. First, 
note that the fermion d.o.s. is constant at the Fermi 
energy, i.e. 



lira Nt(u;,g)=Nf(e F ,x) + — . (B4) 



Next consider the boson d.o.s. 



N b (uj) = d 2 k6 



2U 



ifeCfe 



(B5) 



Clearly N b (uj) = for |w| < A and for w - A « A 
we can expand the dispersion at the bottom of the band 
e k ps |/c 2 to obtain 



JjmJV 6 ( W ) = ±*( W -A)+o(£). (B6) 



Plugging both these expression into Eq. (|B3p , we obtain 

jV(w) = — 6>(w - A) / dilNf (w - A - 0) (B7) 

Nf(e F ,g) 



(u - A)9(uj - A)- 



2/ 



(B8) 
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